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In the presence of cosmic chiral asymmetry, chiral-vorticity and chiral-magnetic effects can play 
an important role in the generation and evolution of magnetic fields in the early universe. We include 
these chiral effects in the magnetic field equations and find solutions under simplifying assumptions. 
Our numerical and analytical results show the presence of an attractor solution in which chiral effects 
produce a strong, narrow, Gaussian peak in the magnetic spectrum and the magnetic field becomes 
maximally helical. The peak in the spectrum shifts to longer length scales and becomes sharper with 
evolution. We also find that the dynamics may become non-linear for certain parameters, pointing 
to the necessity of a more complete analysis. 

H ~> I. INTRODUCTION 

Cosmic relics, such as topological defects and magnetic fields, can be used as probes of the very early universe and 
extremely high energy particle physics. The formation, evolution, and observational signatures of topological defects 
have been studied for over three decades, and constraints have been obtained that have guided particle physics model 
building. The formation of magnetic fields during cosmological phase transitions has also been investigated [1-3]. 
Current observational constraints are relatively weak, allowing for intergalactic magnetic field strengths at the nano 
i-£h Gauss level (see for example, [4-8]). 

An attractive scenario links the generation of magnetic fields to the generation of the observed matter-antimatter 
asymmetry [9-12]. The produced magnetic field carries magnetic helicity that is directly proportional to the baryon 
number density [9, 10] 

C3 h=-K^, (1) 
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y—i where a is the fine structure constant, k « 0.01 [11, 12], and 
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O^N h = - / ofx A • B, (2) 

V Jv W 

is the magnetic helicity density and rib is the baryon number density. The minus sign in the relation (1) is a direct 
cosmological manifestation of CP violation in particle physics that also gives preference to matter over antimatter 
in the universe. The injection of helical magnetic fields into the plasma can transfer magnetic field power to larger 
^vq length scales by the "inverse cascade", providing hope that even small scale magnetic fields from phase transitions 
can grow to astrophysically relevant scales at more recent epochs. Helical magnetic fields can possibly be detected 
r* through various cosmological observations [13-15]. 

Two important chiral effects, called the chiral-vorticity (xw) [16] and chiral-magnetic (x-B) [17] effects, can also 
play a role in the early universe [18] and in QCD [19, 20]. To understand these effects, consider Fig. 1 where we 
show the effect of a magnetic field on, for example, electrons. The magnetic field couples to the magnetic moments 
and tends to align them. Depending on the electric charge of the carrier, the spins are either aligned or anti-aligned 
with the magnetic moment. The helicity eigenstates of the fermions then determine the direction of the momentum 
of the particles, which in turn gives the direction of the electric current due to each species. Taking the left- and 
right-handed electrons and positrons as the four fermion states, which we assume to be massless, we see that the net 
electric current is 



J xB oc [n(e L ) - n(e+)} - [n(e R ) - n(e+)], (3) 

where n(e~[) denotes the number density of e£ and similarly for the other particle species. In terms of the chemical 
potentials for left- and right-handed electrons, the differences within the square brackets are given by /il and hr 
respectively. So J xu cx A/i = fi^ — /i^. The calculation in Ref. [17] gives 
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FIG. 1: Understanding the \B effect. An external magnetic field tends to align the magnetic moments of the four electron 

states - left-right handedness for electron and positron, denoted in the figure as L+, L— , R+, R which implies the shown 

directionalities of the spin, momenta, and electric current due to each state . If the four states are present in unequal numbers, 
net electric current may be induced. 
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FIG. 2: Understanding the x u effect. Vortical fluid flow tends to align the spins of the four electron states which implies the 
shown directionalities of the momenta and electric current due to each state . If the four states are present in unequal numbers, 
net electric current may be induced. 

Similarly, in Fig. 2, we explain the effect, which occurs if the ambient fluid flow has vorticity (u>). Spin-orbit 
coupling tends to align the spins of the fermions; particle helicity then aligns the left-handed states but anti-aligns 
the right-handed states, which leads to the electric currents as shown. Thus, in equilibrium, 

J XU1 oc [n{el) + n(e+)] - [n(e R ) + n(e+)}. (5) 

The presence of non-zero fi^ means that n(e^) ^ n(et) and of that n(e R ) ^ n(ej). However, if /i£ = /ir then 
n(e£) = n(e^) and n(e^) = n(e^), and J xu vanishes. Also if /i£ = — /ir then n(e£) = n(e\) and n(e R ) = n(e^), and 
again J XUJ = 0. So for J xu to be non-vanishing, we need A/i 2 = /i| — n R ^ 0. The exact calculation in Ref. [16] gives 

J*- = ^2 A ^ 2 " . (6) 

where uj = V x v is the fluid vorticity. 

The above expression for J XbJ holds when the left- and right-handed particles and antiparticles are in thermal 
equilibrium at the same temperature. If some of the species are at different temperatures there is an additional 
contribution per species to J XU1 proportional to eT 2 u) where T is the temperature of the particular species [16]. We 
will not consider this situation in the present paper, though it may be important for the contribution of left- and 
right-handed particles, especially neutrinos, to the hypercharge current in the epoch before electroweak symmetry 
breaking. 

The \B and x 00 effects can only lead to a non-zero electric current if there is a disbalance between left- and right- 
handed particles, that is, A/i 7^ 0. Such a disbalance can arise in the early universe from out-of-equilibrium P-violating 
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decays of massive particles. For example, it could arise due to inflaton decay at the end of inflation. The resulting 
values of Ap are not generally suppressed by any small couplings, so large values like A/i/T ~ 0.1 — 1 can easily be 
achieved [21, 22]. A right-left asymmetry is also likely to be created during leptogenesis, if it occurs at some energy 
scale much higher than the electroweak scale. In particular, an asymmetry between left- and right-handed electrons 
will be created because they carry different charges under the electroweak symmetry group SU(2)l X U(1)y- The 
chirality-flipping processes due to nonzero electron mass are suppressed at high temperatures, and this asymmetry is 
preserved until the temperature drops below Tf ~ 80 TeV [23] . If vorticity develops at any temperature higher than 
Tf, the xw effect will operate and create a magnetic field. Since the electroweak symmetry is unbroken, this magnetic 
field will be the U{X) hypercharge magnetic field, and will get converted to electromagnetic magnetic field after the 
electroweak phase transition. 

It is worth noting that chiral effects induce kinetic and magnetic helicities even if these are not present initially. 
For the effect, Fig. 2 shows that momenta tend to be aligned along u), which means that v ■ ui ^ and so kinetic 
helicity is induced. In the case of the x-S effect, if the magnetic field is helical, then the induced current, hence velocity, 
is along the magnetic field and the fluid flow carries kinetic helicity. This fact may be of interest if the kinetic helicity 
survives until recombination because then it can produce parity odd temperature-polarization correlations [24]. 

Once a magnetic field is generated, its subsequent evolution may also be influenced by chiral effects. In particular, 
Joyce and Shaposhnikov [22] have shown that x-B currents can induce an exponential growth of the magnetic field on 
sufficiently large length scales. A helical magnetic field, in turn, back-reacts on the evolution of the left-right disbalance 
Ap. This is described by the chiral anomaly equation, which relates changes in Ap to changes in magnetic helicity 
[22], 

d(Ap) c A adh 

Here, Tp is the chirality-flipping rate and ca ~ 1, Q is the fine structure constant, and T is the temperature at time 
t. Boyarsky et al [25] have shown that as a result the left-right disbalance can survive down to T ~ 10 MeV if helical 
magnetic fields are present, for example due to baryogenesis. 

The plasma equations in the radiation-baryon single fluid approximation are 

d tP +^-(pv) = 0, (8) 
4 4 

3<9t(H - ^x(Vx«) = JxB-Vp, (9) 

VxB = J, (10) 
V x E = -d t B, (11) 

where we have used natural units, H = 1 = c. The fluid density is p, the pressure is p. The displacement current has 
been ignored, as is done in magnetohydrodynamics (MHD) when the flow velocities are small compared to the speed 
of light. Since we are mainly focusing on chiral effects, we have simplified the equations by ignoring the injection 
of magnetic fields by external sources such as sphalerons [10], and the dissipative effects of viscosity. Cosmological 
expansion can be included in the Maxwell equations by going to conformal coordinates, as we will do in Sec. II. 
The electric current is given by the sum of the Ohmic and chiral components 

J = Johm + J X u + JxB- (12) 

The Ohmic component is given by 

Johm = a(E + v X B). (13) 

where a is the electrical conductivity of the plasma. The chiral components, J x b and J X u, are given in Eqs. (4) and 
(6), together with the chiral anomaly equation in the form of Eq. (7). 

In Eq. (11) we can replace E by the currents and with a little algebra we find 

d t B = V x (v x B) + r )£,V 2 B + j u V x u> + lB V x B, (14) 

where 

1 eAp 2 e 2 Ap 

Id = - , = -j-j- , IB = tt^- ( 15 ) 

The terms on the right-hand side of (14) will be referred to as the advection, diffusion, xw, and %-B terms respectively. 
The above set of equations exhibits all the complexities of MHD plus those arising from the chiral anomaly. 
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It is helpful to compare order of magnitudes of the various terms on the right-hand side of Eq. (14). If L denotes a 
length scale of interest, the advection, diffusion, %w and X-B terms are estimated as 

„ , „s vB B „ eAu 2 v „ e 2 AuB , . 

V x (« x B ) ~ — , 7£) V 2 .B ~ , 7,Vxu £=- , 7s V x £? . 16 

Assuming that jUj, ~ //r ~ T and using cr ~ T/e 2 [26, 27], the estimates for the four terms become vB/L, e 2 B/(TL 2 ), 
e 3 vT/L 2 , and e 4 B/L respectively. Now for small flow velocities, v <§; max{e 2 /(TL), e 4 }, or small magnetic fields, 
B <C e T/L, the advection term is subdominant and can be ignored. The second case we consider is when the length 
scale is large and the magnetic field is small: L 3> min{e 2 / (vT) , l/(e 2 T)}, B <C evT 2 . Then the diffusion term can 
be neglected. 

We solve the evolution equations in an expanding spacetime, first numerically in Sec. II C and then analytically in 
Sec. II D. Throughout this analysis we adopt some simplifications. First of all, we do not attempt to solve the fluid 
dynamics equations and represent the velocity field by a mode distribution with a Kolmogorov spectrum and fixed 
phases. (A more realistic representation of turbulent dynamics would include time variation of the phases.) We also 
assume negligible advection; the consistency of the latter assumption is discussed in Sec. II D. Our conclusions are 
summarized and discussed in Sec. III. In particular, we argue that under certain conditions the magnitude and the 
spectrum of the magnetic fields produced by our mechanism are not sensitive to the details of the turbulent velocity 
flow. Our estimates should then be valid, despite the oversimplified treatment of turbulent dynamics. In Appendix A 
we also consider steady-state solutions, i.e. with d t B = 0, in the regimes of negligible advection and of negligible 
diffusion in Minkowski spacetime. 

II. MAGNETIC FIELD GENERATION AND EVOLUTION 

As discussed at the end of Sec. I, there is a range of conditions within which the advection term may be dropped. 
Without the advection term, Eq. (14) is linear in v and B, and a mode expansion converts the equation into a set 
of ordinary differential equations. However, the dynamics is still highly non-trivial because of the mode coupling that 
arises due to the evolution of A/i, as the magnetic helicity on the right-hand side of Eq. (7) is an integral over all 
modes. Here we shall assume that advection is negligible. The validity of this assumption will be discussed at the end 
of Sec. II D. 

A complete analysis, which is beyond the scope of the present paper, would need to include the evolution of the 
velocity field which is governed by the Navier-Stokes equation. Here we will consider the simpler case of incompressible, 
turbulent flows, with a specified distribution of velocities. These assumptions are valid when the magnetic field energy 
density is much smaller than the kinetic energy of the fluid and we can ignore the backreaction of the magnetic field 
on the fluid flow. 

In a cosmological setting, we assume a flat Robertson- Walker metric 

ds 2 = R 2 {r l ){-drj 2 + tfydxW), (17) 

where R(rj) is the scale factor. It is convenient to choose R(r]) to have dimensions of length, and rj, x l to be di- 
mensionless. The scale factor is related to cosmic temperature by R = 1/T. In the radiation dominated epoch, we 
also define the conformal time as r\ — M*/T, where M* = y / 90/87r^-Mp and g„ ~ 100 is the effective number of 
relativistic degrees of freedom. With this normalization, we find r\ ~ 0.1 when T <~ Mp and, for example, r\ ~ 10 27 at 
matter-radiation equality when T ~ 1 eV. We also define comoving variables 

B c = R 2 (v)B(v), A/V = R(fj)An, (18) 

and the comoving electrical conductivity is given by [27], 

cr c = R{rj)o = 70. (19) 

In these comoving variables, the magnetic field evolution Eq. (14) takes the form 

d v B c = V c x (v c x B c ) + lD( XJ 2 c B c + 7wc V c x (V c x v c ) + 7Sc V c x B Cl (20) 

where 7 £> c = l/u c , 7aJC = eA/i 2 /47r 2 cr c , jbc = e 2 A/^ c /27r 2 cr c , V c is differentiation with respect to the spatial metric 
5ij and v c is the comoving velocity. Note that v c — v by definition. 

From Eq. (7), the evolution of the comoving chemical potential is given by 

d ^ lc = -c A ad n h c - T F A^ C , (21) 
ai] 

where h c is R?h and Tp is the chirality-flipping rate that we can ignore in the early universe T > 80 TeV [23]. In 
what follows we will always work with comoving quantities and will omit the subscript c for convenience. 
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A. Decomposition into modes 

In order to solve for the evolution of the magnetic field with a chemical potential, we decompose the vector fields 
in the modes, Q i , which are divergence- free eigenf unctions of the Laplacian operator in comoving coordinates, 

Q±(fc)= ei(fc)±fe 2 (fc) exp( , fca;)) (22) 
v2 

where e% ~ k/k and (ei, 62,63) form a right-handed, orthonormal triad of unit vectors. Then, V • Q ± = and 
V x Q ± = ±kQ ± , and we also take Q ± *(-fc) = Q ± (+k). 

The velocity field of the incompressible (V • v = 0) fluid is now decomposed in modes 

/d 3 k 
—^3 [v + (r,,k)Q+(k)+v-(r,,k)Q-(k)] . (23) 

The fluid kinetic energy density is given by 



(\v( V) x)\ 2 ) = ^ JdlogkE v (r 1 ,k)=p r J^[\v + (r ] ,k)\ 2 + \v-(r ] ,k)\ 2 ], (24) 



where p r is the radiation density, and we have taken statistically isotropic correlators, 

(S^fa.fc^fap)) - l^fo.fc)! 2 (2^) 3 <5 (3) (fc-p) (25) 
(i) + *(?7,fe){;-(?7,p)) = (v-*(r},k)v + [r),p)) = 0. (26) 



Proceeding in the same way, the fluid helicity is 

(v.u) = J 'dlogk H v (r,,k) = J^[\v + (v,k)\ 2 -\v-(v,k)\ 2 ]. (27) 



The magnetic field is similarly decomposed 

r rl 3 k r 1 
B( V ,x) = J ^ [ J B+(7 ? ,fc)g+(fe) + S-( ?7 ,fc)Q-(fc)J . (28) 

where the scalar mode is absent because V • B — 0. The ensemble average of the magnetic field energy density is 

1 r r k 2 dk 

2 <|B(,7,x)| 2 > = J dlogk E B ( V ,k) = J j^ 2 [\B + { V M 2 + \B-{vM 2 ]- (29) 

where, just as for the velocity field, 

(fl** fa, P)) = l^fofc)! 2 (2Tr) 3 5 (3 \k-p) (30) 
{B+*{ v ,k)B-{r 1 ,p)) = (B-*(r,,k)B+(r l , P ))=0. (31) 

The magnetic field helicity density is 

(A-B) = JdlogkH B (r,,k) = |^[|B + (»?,fc)| 2 -|B-fa,*)| 2 ]. (32) 
and from Eq. (2) we have 

(fc> = <A.B>. (33) 

The MHD equation, Eq. (20), without the advection term, can be decomposed into equations for the modes -B 

d ri B+ = {- lD k 2 + lB k)B+ + lul k 2 i+, (34) 
d v B~ = (-j D k 2 -j B k)B- + lu k 2 v~. (35) 

We multiply the first equation by B + * and the second by B~* and take ensemble averages to get 

d v \B+\ 2 = 2(- lD k 2 + lB k)\B+\ 2 +2 lu k 2 (B+*i+), (36) 
d v \B-\ 2 = 2{- lD k 2 - lB k)\B-\ 2 +2 lu k 2 {B-*v~). (37) 
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The magnetic field is zero initially and only the %w term is important on the right-hand side. Hence at early times, 
the solution to Eqs. (34) and (35) is 

B ± ( v ,k)= lu k 2 f dr,' e ± (r 7 / ,fe) (38) 

Jvo 

Note that 7 W is, in principle, a function of time since A/x can vary due to the chiral anomaly relation. However, the 
Xw term cannot change the helicity of the magnetic field as can be checked by taking the difference of (36) and (37), 
and so 7^ can be assumed constant and taken out of the integral. Then for the %w term we get 

(B ± *(r 1) k)v ± (r 1) k'))=j u; k 2 f drf (»/,*)«*(»/, *')> (39) 

and we need the unequal time correlator for the velocity field. 

We expect the fluid velocity at any k mode to be correlated on the eddy turnover time scale 27r/fcw(ry, k), where 
v(ji, k) is the fluid velocity, and to be uncorrelated on longer time scales. This suggests 

(v ± *( v / ,k)i ± ( v ,k / )) = (i ± ( v ,k) 2 ) (2nf6^(k-k'), for \v-rf\ < 1 ~^r (40) 

and 

2-7T 

(v ± *( V ',k)v ± (r ] ,k')) = 0, for |r7 - t/| > — — - (41) 

where v(t], k) is the fluid velocity at length scale 2-K/k at time 77. (We will relate v(t], k) to v (f), k) in Eq. (55) below.) 
Inserting these unequal time correlators in Eq. (39) gives 

(S^fa, *0> = 7u k 2 f( V ,k) |«±| 2 (27r) 3 S^(k-k / ) . (42) 

where f(j), k) = r\ — r; for r/ — t/q < 2ir/(kv), and f(r], k) = 2ir/(kv) for 77 — 770 > 2ir/(kv). An example of a smooth 
function with the same asymptotic properties is 

/( r? ,fc)=^tanh^(r ? - 77o )) (43) 

where S ~ 1 is a fudge factor. We will adopt this choice for /(tj, k) in our numerical calculations. 

The form of the correlator (42) is based on (39) which is valid at early times when the \w effect dominates. The 
correlator will not be valid at later times when the %u; effect becomes sub-dominant, but then the form of the correlator 
is also not important. Hence we can adopt the form (42) for all times. 

So the MHD equations become 

9,|S+| 2 - 2(- lD k 2 + lB k)\B+\ 2 + 2^k i f( V ,k)\v+\ 2 (44) 
d.lB-f = 2(- lD k 2 - lB k)\B-\ 2 + 2 1 2 u k i f( V} k)\v-\ 2 (45) 

These equations are linear in the quadratic variables l-B*! 2 . 

From Eq. (21), the evolution of the ensemble averaged comoving chemical potential is given by 

lif = ~ CA( * I 2^" [\B + \ 2 -\B-\ 2 ]-r F A^ . (46) 
Eqs. (44), (45) and (46) are the dynamical equations that we will need to solve after specifying the velocity modes, 



B. Fluid Velocity 

We consider the scenario of a strong first order phase transition at an early epoch, much earlier than the electroweak 
phase transition. As bubbles of the new phase grow and collide, the cosmological medium gets pushed and turbulence 
is generated. The radius of the bubbles at completion of the phase transition defines the length scale at which the 
fluid is being driven. Let us call this length scale Cb- In a non-expanding spacetime, on length scales smaller than ^ we 
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expect a Kolmogorov distribution of velocities, and on larger scales, if we assume a random superposition of constant 
velocity domains of size we expect a white noise spectrum. In a cosmological setting, we also need to account for 
Hubble expansion. The largest scale relevant for turbulent flows will be the size of the eddy that circulates at least 
once per Hubble time. This length scale is Ch ~ vrj where v is the flow velocity driven by the bubbles. Then the 
inertial scale is set by the smaller of and (h and will be denoted by (j. 
The fluid velocity spectrum is given by a power law 



Ev ^= v ^{m) ' (47) 

where fej denotes the inertial wavenumber, is the fluid velocity at this scale. On small length scales, k > hi, the 
Kolmogorov spectrum gives n = —2/3, whereas we take a white noise spectrum, n = 3, on large length scales, k < k t . 
(The assumption of white noise on large length scales is not crucial for the main qualitative features in the evolution 
as these all occur in the Kolmogorov part of the spectrum.) From Eq. (47) we get the velocity at wavenumber k, 

( k \ n/2 

<V,k)^M{ m ) ■ (48) 

In the radiation era, and on large length scales with k < fcj(^), the velocity field does not change (e.g. Sec. 7.3.2 of 
[28]). So 

v(r],k) = v(r)o,k), k < fctfa), (49) 
where 770 is the time of the phase transition. Combining Eqs. (48) and (49) with n — 3, we obtain 

vM=v i ( Vo )(P^) 32 ■ (50) 



Now the inertial wavenumber, fcj, is defined by 



When inserted in Eq. (50) we get 



Plugging back into (51) gives [29] 



h(v) = (51) 
u< (77)77 



, 3/5 
770 x 



vM=Vi(vo) Hr ■ ( 52 ) 



2/5 



h(v) = HVo) If) ■ (53) 



With these expressions the fluid velocity in Eq. (48) can be written as 



x (3-n)/5 / , x nil 

770 V / k x 



We will assume that there is no kinetic helicity in the velocity flow and from Eq. (24) write 

|u+(?7,fc)| = \v-{y\,k)\ = i:k- i/2 v{ri,k). (55) 



C. Numerical Evolution 



We will now solve Eqs. (44), (45) and (46) with the velocity spectra as given in Eqs. (55). We will consider a 
few different values of 7^(770) and take 770 = M*/Tq where T = 10 10 GeV is the temperature at which the phase 
transition occurs and = 6.6 x 10 17 GeV. At the high temperatures we are considering, e should be the Abelian 
(hypercharge) gauge coupling constant but, for numerical purposes, we take it to be given by e = \[^Ka in 7^ and 
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FIG. 3: Plots of the magnetic energy spectrum at T = 10 9 GeV for different values of the peak velocity and fudge factor S (see 
Eq. (43)) exhibit very similar features. 
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FIG. 4: Evolution when the initial velocity flow is non-helical with Vi(r]o) = 10 -4 and the initial chemical potential is A/io = 1. 
The left panel shows the power spectra, Eb and Ev, which are normalized to the comoving radiation energy density, <;,7r 3 /15. 
The thick lines are for Eb and the thin lines are for Ev- The spectra are shown at the three different times, T = 10 9 GeV 
(dashed), T = 10 7 GeV (dotted) and T = 10 5 GeV (solid). The thin dotted-dashed line is the initial velocity power spectrum at 
T = 10 10 GeV while the initial magnetic field vanishes. Similarly, the right panel shows the evolution of the magnetic helicity 
spectrum. 



7b, where a = 1/137. The comoving electrical conductivity is a = 70 as in Eq. (19). We consider several different 
values for the chemical potential at t]q but focus on AfiQ — 1.0 = A/!q. The initial magnetic field is taken to vanish in 
all cases, as is the flipping rate Tp since this is small at temperatures above ~ 80 TeV [23]. Note that the equations 
of motion are independent of the initial epoch of turbulence but the velocity field in Eq. (54) explicitly contains 770, 
and hence T . 

We start by showing the magnetic energy spectrum at a fixed time (T = 10 9 GeV) for several different values of 
the peak velocity and parameter S (defined in (43)) in Fig. 3. The plots show that a sharp peak in the spectrum 
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FIG. 5: Evolution of the total helicity (left panel) and the chemical potential (right panel), when the initial velocity flow is 
non-helical and the initial magnetic field vanishes, for several different parameters. 
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FIG. 6: Evolution of magnetic energy density with temperature for three different values of and with Vi(rjo) — 10 



develops and its position and shape are not very sensitive to the input parameters. This can be understood from the 
evolution equations (44) and (45). For jb > 0, we find that B~ stays small and only B + contributes to the magnetic 
energy density and helicity. At early times, the velocity field acts as a source term for B + and the first term on the 
right-hand side of Eq. (44) is negligible. However, with evolution, the first term becomes more important and the 
coefficient is negative for k > ^b/id and positive for k < ^b/id- This change in sign implies a peak in the spectrum 
at k Jb/jd ~ e 2 A/x/27r 2 , or on a comoving length scale 47r 3 /(e 2 A/i). The position of the peak at k w Jb/jd also 
agrees with our numerical results. We will give a more precise analytical derivation of the peak location in Sec. II D. 

The left panel of Fig. 4 represents the evolution of Es(k) and Ey(k), and the right panel shows the evolution of 
Hs{k), for A/io = 1 and Vi = 10~ 4 . As is evident in Fig. 4, the peak position shifts toward large length scales with 
evolution. This can be understood by noting that the peak position oc A/i, and A/i is a decreasing function of time. 

The evolution of the chemical potential and magnetic helicity are shown in the two panels in Fig. 5. The difference 
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of chemical potential, A/x, decays and becomes smaller, while the total helicity, h, evolves to its asymptotic value 

h* = A/xo/acA- (56) 

The main contribution to the total magnetic helicity comes from the peak of the helicity spectrum. At this peak, 
the diffusion and X-B terms in Eq. (44) cancel each other. Further evolution occurs since the chemical potential is 
changing and this shifts the peak position. 

The chemical potential difference A/x has an "attractor" solution - different initial values, A/xq, all evolve to decay 
along a common trajectory with T 1 / 2 slope (see Fig. 5). The attractor behavior we observe is stronger than the 
"tracking" observed in [25] where the authors find constant helicity at late times in the case of a single magnetic field 
mode. Joyce and Shaposhnikov [22] also observed the T 1 / 2 decay of A/x; in addition, we find that the coefficient of the 
decay is independent of the initial value of A/x. Eventually A/x becomes very small and the magnetic helicity reaches 
its asymptotic values (56) as given by the chiral anomaly equation. We will obtain a better understanding of these 
features in the next section. 

In Fig. 6 we plot the magnetic energy density as a function of temperature. During the attractor evolution, this 
gives pb oc T 1 / 2 . 

Up to now, we have neglected the advection term in the magnetic field evolution. This assumption is certainly valid 
in the early stages of the evolution when the magnetic field strength is very small and the chemical potential is large. 
At later times, however, we might expect the advection term to become important. To clarify this issue we compare 
the advection and terms at the peak position after we have analytically understood the evolution in Sec. II D. 

Notice in Fig. 4 that the magnetic energy density at the peak is larger than the kinetic energy density. One might 
expect backreaction of the magnetic fields on the fluid flow to become important in this situation, and the magnetic 
and kinetic spectra to asymptotically approach equipartition. However, in our case, except for a brief initial phase, the 
current is proportional to the magnetic field and so the Lorentz force, j x B, vanishes. Thus there is no back-reaction 
of the magnetic field on the fluid flow. 



D. Analytical Understanding 

In this section we will be able to understand most of the features of the numerical evolution. Our analysis also 
applies to the scenarios considered in Refs. [22, 25] since the only difference is in the early stages of the evolution. 
In Ref. [22], A/xo 7^ and a non-zero initial magnetic field are assumed. In Ref. [25], A/xo = and the initial helical 
seed field is injected during baryogenesis. In the present work A/x ^ and a non-zero seed field is generated due to 
turbulence and the x^ effect. 

The equations to be evolved are (44), (45) and (46). Initially, B ± = and so only the x w terms are important. 
Then 

1^(77, k)\ 2 « 2fc 4 r di) ltf(v,k) 1^(77, k)\ 2 (57) 
J no 

Note that 7^ has not been pulled out of the time integration because it is proportional to A/x and this must be evolved 
according to (46). However, at these early stages, B + « B~ and the right-hand side of (46) vanishes, and so A/x and 
7^ are constant. This is also seen in Fig. 5 where the curves of A/x are flat at the highest temperatures. 

The x w dominated evolution continues until the \B term starts to become important. The transition occurs at 
time 771 such that 

rm 

lB k ^/(r/,fc)|« ± (7 7 ,fc)| 2 « f{rn,k)\»HmM 2 (58) 

Jr/o 

where v ± are given in Eq. (55). For 77 > 771 , we can ignore the x^ term. 

In the Kolmogorov part of the spectrum, Eq. (54) gives v ± oc 7y -11 / 15 and, assuming 2ttjb < 1, the integral on the 
left-hand side of Eq. (58) is dominated by the contribution from 77 ~ 2n/k. A simple estimate then gives 

1 2tt 

**( 27r7fl )lB/22T (59) 

and this introduces an extra k dependence in B ± (j]i,k). However, we find B ± (rji(k),k) oc fc -1 / 10 in the Kolmogorov 
part of the spectrum, and this weak dependence will be ignored in what follows. Also, as in Fig. 4, the evolution 
will lead to a sharp peak in the spectrum, which is the most interesting and dominant feature in the spectrum. In 



11 



discussing this peak, we can treat 771 as being independent of k, since the Kolmogorov spectrum is relatively flat over 
the width of the peak. So we will simply set 771 = from now on. This simplification will not affect the gross features 
of the evolution. 

Ignoring the x w term, the equation for B ± can be written as 

= -7^(^^)15*1 (60) 

with 

7s e 2 A/x 

k P = — = -t^t ■ ( 61 ) 

Eq. (60) can be solved and the solution written as the product of an initial amplitude, an exponential amplification 
factor common to all modes, and a Gaussian spectrum, 

^±(77, fc)| = \B±\ cxp( lD K 2 r,) exp [- 7D (k T K p f 77] (62) 

where the peak of the Gaussian for B + is at 



1 



n 



k p = I dr ) k p ( 63 ) 

The restriction k > implies that B~ is given by the tail of a Gaussian which is peaked at k = —K p . Therefore B~ 
can be ignored except for 77 pa 0. 

Thus the spectrum of B + is a Gaussian of width 

Afc Pa — *— (64) 

which is consistent with the width of the peaks in Fig. 4. It is also interesting to note that the peak gets narrower 
with time. 

Now we can estimate the magnetic helicity 

^- 2 k(\B^-\B-\ 2 ) (65) 
PS —K v e 2 ^ K > (66) 

where A = |-B^~| 2 / \f2-K 2 ^£,. The magnetic energy density is obtained from Eq. (29) 

1 f dk 

-(\B(r h x)\ 2 ) = J — fc 2 [(|B+(77,fc)| 2 ) + (|S-(77,fc)| 2 )] 

* ^, K > 2lDK> (6?) 
From the chiral anomaly equation, we also have 

A/x = A/io — CAah (68) 

« A M0 - °-^K p e^> (69) 

As discussed above, during the xcj dominated phase, the helicity stays near zero and A/x pa A/i . Then fc p pa 
e 2 A/io/27r 2 and -K~ p pa fc p /2 pa constant. In a time pa [7£)(A/xo) 2 ] _1 ; the exponential in (69) starts to become important. 
Then the magnetic helicity increases exponentially fast and A/j rapidly decreases. 

The exponential growth of the helicity terminates when it gets close to the asymptotic value /x*, given by Eq. (56). 
The subsequent evolution can be investigated by setting h pa in Eq. (66). In order for h to saturate at h*, the 
exponent in Eq. (66) has to become nearly constant, K 2 ij pa const, or 

K p cx 77" 1/2 . (70) 
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Then Eqs. (63) and (61) imply 

k p « K p (71) 

and 

A/i cx ?r 1/2 oc T 1 / 2 (72) 

which is the attractor solution seen numerically. 
To make these estimates more quantitative, we set 

K p ^k p = ^^CT^, (73) 

where C is a constant (or, more exactly, a slowly varying function) to be determined. Substituting this into (66) and 
using T = M t /r], we obtain 

C^{2 lD AU)- 1 ' 2 {\nX) 1 / 2 , (74) 

where 

X = {h*/AK p )y/rj oc 77. (75) 

This shows that C has a weak, logarithmic dependence on 77. 

Comparing Eqs. (66) and (67) and using (73), we can estimate the magnetic energy density in the attractor regime, 

i<|B(^x)| a >«i^«^rV». (76) 

Note that, apart from the logarithmic factor, the coefficient C in (73) is a constant independent of A//o and of the 
turbulent velocity spectrum v (r),k). Thus Eq. (76) shows that the attractor magnetic energy density is (roughly) 
proportional to the initial chiral disbalance A^t and is largely independent of v (r], k). The characteristic length scale 
of the magnetic field, 

K- 1 a C- l T- x l 2 , (77) 

is not sensitive to any of the input parameters. 

To get numerical estimates, we evaluated X and C at the lower end of our temperature range, T — 10 5 GeV. This 
is done by solving the transcendental equation (74), after substituting expressions for A (below (66)), B+ (from (57)), 
771 (from (58)), and other relevant quantities. For A/jq = 1-0, this gives X rj 3.8 x 10 16 , C rj 4.5 x 10 -8 , in reasonable 
agreement with our numerical results. We have verified that for the range of T and A/io that we considered here C 
does not vary from this value by more than 0.8%. 

We now compare the advection term to the \B term at the peak position. The estimate in Eq. (16) shows that the 
advection term is sub-dominant if 

v x < 7s (78) 
where v x denotes the velocity in physical space. We use Eq. (24) to go to momentum space 

v x ~v(r],k) (79) 
and v(r), k) is given in Eq. (54). Now we set k = K p = CT 1 / 2 , r] = M*/T in (54) to get 

v{r, ' Kp)=Vio {^) UJ (80) 

where = fj(?7o) and n = —2/3. Inserting numbers in GeV units Vio = 10 -4 , C = 4.5 x 10 -8 , M* = 6.6 x 10 17 , 
T = 10 10 we get 

rp \ 17/30 



v(rj,K p ) = 6xl0- a [-) (81) 



13 



We also have 

The temperature dependence of v and 7^ are very similar and so whatever relation holds at the initial time, also 
holds for all (relevant) later times. For a range of parameters we have used in the numerical analysis, v < 73, and 
so the advection term is only marginally important. (From Eq. (74) we see that C is only sensitive to the parameter 
7d = 1/cr.) A smaller value of Ujo, or a larger value of 75 due to several chiral particle species contributing to the \B 
effect, or a higher initial temperature T , can further ensure that the advection term is unimportant. We also note 
that the advection term becomes negligible when turbulence is eventually dissipated on the relevant scales. However, 
if the condition v <C 75 is not satisfied, the advection term will be important and a full numerical analysis will be 
necessary 

To summarize our understanding, at very early times, only the x^° effect is important and this generates magnetic 
fields from the assumed turbulence. As these magnetic fields get stronger, the dissipation and \B effects become more 
important than the %w effect. These two terms take the initial spectrum produced by the xw effect, and introduce 
a Gaussian peak in the spectrum at k = K p whose amplitude grows exponentially and with width that decreases as 
1/y/r}. At the same time the magnetic helicity h increases exponentially, the magnetic field becomes maximally helical, 
and A/i rapidly decreases. The rapid growth of h terminates as it approaches the asymptotic value ft.*, and the system 
enters the attractor regime, in which K p w k p oc T 1 / 2 and Afj, oc T 1 / 2 . This means that the peak of the spectrum 
continues to evolve to larger length scales. In the attractor evolution, the magnetic helicity stays approximately 
constant while the magnetic energy density falls off as T 1 / 2 . 



III. CONCLUSIONS 



Maximal parity violation in the standard electroweak model suggests that there may have been an asymmetry in the 
distribution of left- and right-handed chiral fermions in the early universe, so that the cosmic medium was a "chiral 
plasma" . Such a chiral asymmetry has implications for the generation of magnetic fields via the chiral-magnetic and 
chiral- vorticity effects [16, 17]. The latter effect leads to the production of magnetic fields during a turbulent phase in 
the chiral plasma, while the chiral-magnetic effect leads to amplification of the magnetic field [22, 25]. These effects 
are in addition to the generation of magnetic fields during phase transitions [1, 2] and during sphaleron processes 
responsible for baryogenesis [2, 3, 9-12]. 

In this paper we have focused on the consequences of the xu> and x^ effects. We have investigated the time- 
dependent solution numerically in an expanding universe but with two assumptions. First we have considered negligible 
backreaction of the magnetic field on the fluid velocity, and second we have assumed that the advection term is 
unimportant. These assumptions simplify the analysis because they linearize the MHD equation. However, the analysis 
is still highly non-trivial because the chiral anomaly equation connects the evolution of the chemical potentials to the 
total magnetic helicity which is a sum over all magnetic modes. 

Neglect of back-reaction is justified by the fact that maximally helical magnetic fields generated by our mechanism 
are nearly force-free. We have also verified that it is possible to choose parameters such that the advection term 
is unimportant for the evolution of the magnetic field at the peak of the spectrum. We leave the inclusion of the 
advection term, and perhaps a full-blown magneto-hydrodynamic analysis, for future work. 

Our numerical solutions show that a magnetic field is generated due to the \io effect in a turbulent chiral plasma. 
As the field gets stronger, the %B effect becomes important and leads to rapid amplification of the field on a preferred 
length scale given by Ib ~ (KpT)^ 1 ~ 47r 3 /(e 2 A/i). 1 The x^ effect becomes insignificant at this stage. The amplifica- 
tion period ends when dissipation (magnetic diffusion) becomes important, and the evolution approaches an attractor 
regime in which the chemical potential difference decreases as T 1 / 2 cx i™ 1 / 4 , with a corresponding shift of power in 
the magnetic field to larger length scales. The magnetic helicity density h remains nearly constant in this regime, and 
the magnetic energy density ps decreases as T 1 / 2 . The attractor solution for psit) is proportional to the initial chiral 
disbalance Ap /T Ql where T is the temperature at the phase transition, and has only a weak logarithmic dependence 
on the magnitude of velocities and on the spectrum of turbulent motion of the plasma. The coherence length scale of 
the field is largely independent of both A^q/Tq and the turbulence spectrum. 

The attractor regime ends at Tp ~ 80 TeV, when helicity flipping becomes important. The magnetic energy density 
at this point is ps ~ 10 _3 (A/i /To) p r where p r is the energy density in radiation, and the coherence length is 



In this section we use the physical, rather than 'comoving' value for the chemical potential difference A/i. 
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Ib ~ 10 5 Tp , which corresponds to the present comoving scale of 10 4 cm. At T < Tp, chirality flipping will reduce 
A/i/T and shift the peak to yet larger length scales as in Ref. [25]. When the temperature drops to ~ 0.1 MeV, 
electron-positron pairs annihilate, resulting in a sudden decrease of the electrical conductivity, cr, by a factor ~ 10 -9 . 
Analysis similar to that in Sec. II D suggests that a sudden decrease in a at a time when Afi = will decrease K p 
and thus increase the coherence scale. In addition, the width of the Gaussian peak in the spectrum will also increase. 
Further growth of the coherence length may be caused by the inverse cascade in the decaying MHD turbulence (see, 
e.g., [30, 31] and references therein). All these effects combined could shift the present coherence length into the range 
of astrophysical interest. We hope to discuss the evolution of the spectrum through the various cosmological epochs 
separately. 

In addition to the time-dependent analysis we have solved the MHD equations for a steady state solution. This 
analysis points to the estimate in Eq. (A15) for the initial field generated by the cosmic turbulence. 
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Appendix A: Steady state solutions 

To obtain the steady state solutions, we set d t B — in Eq. (14) and assume that the velocity field is stationary, 
dtv = 0. 



1. Negligible advection 

We first consider the case when the advection term, V x (v x B), is small compared to the other terms. Then 
Eq. (14) reduces to 

7dV 2 B + 7 u Vxu + 7 B VxB = 0, (Al) 

or 

V x[- lD V x B + 1uj uj + 1b B]=Q, (A2) 

Then 

-7flVxB + 7[J w + 7 B B = V^ (A3) 

where <f> is any scalar function. However, the divergence of the left-hand side vanishes and hence we find 4> = 0. So 
the equation reduces to 



In Fourier space 



71)7x5 + 7^ + 7^ = 0. (A4) 



B = ' Bte ""'' <A5) 



Inserting in Eq. (A4) gives 



- ij D k x B k + i^k x v k + j B B k = 0, (A7) 
The dot product of this equation with k gives 

k B k = 0, (A8) 
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and a dot product with v k gives 

{-ilDVk x k + "/ B v k ) ■ B k = 0. (A9) 

Therefore 

B k = F(k)k x (-i~f D v k x k + ~/ B v k ). (A10) 
To determine the scalar function F(k), insert into Eq. (A7) and obtain 



F = ^°W' (A11) 



Therefore 



B k = — 2 1— k x (~ij D v k xk + ~/ B v k ). (A12) 

1b — 1d k 

If we also assume that the fluid is incompressible, we get 

B k = 2 ~ % \ [ lB k xv k - i lD k 2 v k ] , i£k-v k =0. (A13) 

1b lD k 

We now return to our approximation in this section, namely to ignore the advection term. As discussed at the end 
of the previous section, this is justified for small fluid flow velocities and small length scales. Our solution shows that 
there is another situation when it is fair to ignore the advection term. Suppose the spectrum of velocities is dominated 
by a single mode. Then the advection term in Fourier space is k x (v k x B k ) and this vanishes when we use the 
solution in Eq. (A12). 

The solution (A12) simplifies in the large wavelength limit, k — > 0. Then 

Bk = -i-^k x v k , (A14) 

IB 

or 

B = u> = u>. (A15) 

IB 2e K ' 

This can also be seen directly from Eq. (14) because we have discarded the advection term and the dispersion term 
becomes negligible on large length scales. So the steady state solution is given by 



J X u + J x b = 0, (A16) 



which leads to Eq. (A15). 

To get an estimate for the magnetic field strength with coherence scale L, we assume 



This gives 



Kl ~ V f- (A17) 



B~-^-^v L . (A18) 
ZeL 



In a time-dependent situation with B = initially, we can expect a field of strength (A15) to develop on a time 
scale 

r L ~ L/7 B . (A19) 

This estimate follows by taking the ratio of the x w term in the evolution equation, ~ "f^u/L, and the asymptotic 
value of the magnetic field ~ (7^/75)^. In a turbulent flow, individual eddies are not expected to persist much longer 
than a single revolution time, ~ w^ 1 . This is longer than tl, provided that 

v l <1b- (A20) 
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Note that in this regime, the advection term in (14) is small compared to the term, so the neglect of advection is 
justified. 

If Vjj is given by the Kolmogorov spectrum 



Vl oc v t 



1/3 



(A21) 



where vt is the velocity on a cosmological length scale, and assuming that Vt ~ 1 and (hl + I^r) ~ T, we have 

R T 

2ei 2/3 t l/3 ' 

Inserting numbers suitable for the electroweak scale, we obtain the magnetic field on the length scale L, 

10 18 



(A22) 



Br 



(LT EW ) 2 / 3 



G. 



(A23) 



where T EW ~ 100 GeV. 



2. Negligible diffusion 



We now find the steady-state solution with the advection term but without the diffusion term. Now the steady 
state equation is 



which leads to 



V x (v x B) + 7 W V x u> + 7bV x B = 0, 



v x B + j^lj + 7b-B = V0, 



(A24) 
(A25) 



where cf> can be any function, including a (scalar) function of the flow velocity e.g. v 2 . Unlike the case of diffusion 
domination, the divergence of the left-hand side does not vanish and we cannot argue away <fi. 
By taking scalar and vector products of v with Eq. (A24) we find the solution 



B 



7u 



1b{j 2 b +v 2 ) 



[(v ■ u)')v - -) B v x J + 7|w'] , 



where 



oj = u> Vc 



To determine (f> we impose V • B = 0. This gives the equation for 



1 + u 2 3 



where u — v/jb, & = <fi/jB and 











V x ( 


' u )} 

l + u 2 ) 






(«• 


Sl)u- 


U X ft + 




n 


- u 2 



• V* = 5, 



with $1 = V x u. The parameters 7^ and 75 in the above equations are defined in Eq. (15). 
Eq. (A28) can be re-written in terms of the metric 

g lJ = (l + u 2 )(^'+«V). 



(A26) 



(A27) 



(A28) 



(A29) 



(A30) 



Then the equation looks like a Klein-Gordon equation in a curved background with an O • V$ coupling and a source 
term. 

The solution of (A28) can be inserted into the above expression for B in Eq. (A26) and that will be the steady 
state solution in the case of negligible diffusion. 
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The advection term is important for large velocities, in particular, for t) > 7^, where jb is the parameter in the 
\B effect. In this case, the second and third terms in the square bracket in Eq. (A26) are suppressed by powers of jb 
and we get the order-of-magnitude estimate 

Bk-^v-uj v, (A31) 

IB 

where v denotes a unit vector in the direction of v. The presence of kinetic helicity is crucial for this estimate to be 
non- vanishing. Assuming helical flow with v ■ uj ~ vlo 1 we get 

B~-1^uj, (A32) 
7b 

which coincides with the order of magnitude estimate in the case of negligible advection, Eq. (A15). 
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